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ABSTRACT 

Some radio pulsars show clear 'drifting subpulses', in which subpulses are seen to drift in pulse longitude in a systematic pattern. 
Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the 
subpulse period (P3) is constant on timescales of days, months and years, and between 14 — 5100 MHz. Despite this, the shapes of 
the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through 
the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of ~ 180° in subpulse phase is observed 
at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GET, WSRT and 
Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step 
increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate 
driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does 
at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be 
explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical 
processes happening within the pulsar magnetosphere. 

Key words, pulsars: general - pulsars: individual: PSR B0809-I-74 - telescopes: LOFAR 
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1. Introduction 

PSR B0809+74 is a nearby (0.43 kpc, 
pulsar, whose proximity to Earth and relative brightness 



Brisken et al. 2002) 



quencies between ~ 12 MHz and 


~ 10 GHz (see, e.g. Backer 


et al.|1975 Bruk et al.'1986; Proszynski & Wolszczan 1986| 


Seiradakis et al.||1995. Edwards & Stappers,i2003; .Rosen] 


& Demorest 2011 Gajjar et al. 


2012). The broadband de- 


tectability or the source, as we 


1 as the fact that the pul- 



sar exhi bits many interesting features su ch as 'drifting sub 



'microstructure' ( Craft et al.||1968| ), makes it an excellent 



pulses' dPrake k Craft"l968'), 'nulling' |Backer||1970|) and 

lake 



source for studying the pulsar magnetosphere and the pul- 
sar emission mechanism. Here we focus on PSR B0809-|-74's 
drifting subpulses. 

Drifting subpulses were first identified by [Drake fc Craft] 
( |1968| ), who noticed that the single-pulse components ('sub- 
pulses') of some pulsars are not randomly distributed across 
the emission region. Instead, they appear to drift through 
the pulse window in a well-defined pattern. The drift rate 
of subpulses i s normally characte rised by three parameters 
(as defined by Sutton et al.|1970 ): Pi, the pulsar's rotation 



period; P2, the interval between subpulses; and P3, the sep- 
aration between two driftbands. 

Initially, drifting subpulses were attributed to oscilla- 
tions of the neutron star surface. This ide a was proposed 
with the discovery of drifting sub pulses ([Drake fc Craft 
1968) and was further developed by van Horn ( |1980p , who 
showed that the surface of a neutron star can support 
p-mode oscillations, vibrations caused by pressure waves. 
These surface oscillations modulate the pulsed emission, 
and if the waves are in a beat frequency with the rotation 
period then the oscillations are observed as drifting sub- 
pulses. 

Surface oscillation models were widely abandoned in 
favour of the so-called 'carousel' model (Ruderman & 



Sutherland 1975 ) , because of their failure in explaining 
several features in the observed drifting patterns - these 
i nclude the conser vation of pulse longitude across nulls 
( Unwin et al.|1978), the h igh degree of stability of the drift 
bands ( Staelin et al. '197(7) , and the changing drift-rate after 
a null (Lyne & Ashworth 1983). In the carousel model, ra- 
dio emission is generated from discrete locations ('sparks') 
in the plasma of the pulsar magnetosphere. The configura- 



1 



Hassall et al.: Frequency-dependent Delay from the Pulsar Magnetosphere 



tion of the sparks with respect to each other remains fixed, 
but the whole sparking region rotates around the magnetic 
pole, like a carousel. This rotation (caused by E x B drift) 
moves the sparks slightly between pulses, so the subpulses 
appear to drift. The rate of rotation has been measured for 
a few pulsars, allowing d etailed maps of the polar c ap to 
be made (see for example, Deshpande fc Rankin|200T| ), and 
modified versions of this model have been used to try to 
explain all of the features listed above (see e.g. [Filippenko 
&: Radhakrishnan|1982[[Gil fc Sendyk|2000l [GiTet al.|2003t 



ously unp ublished 4850 MHz dat a from the Effclsberg tele- 
scope (see Lazaridis et al. 2008 for details of the system^ 



2008). 

More recently however, it has been shown that drifting 
subpulses can indeed be explained by surface oscillations, 
as long as the oscillations a re centred on the magnetic axis . 



and not the rotation axis (Clemens & Rosen 2004 20081 



These 'non-radial' oscillations address many of the prob- 
lems mentioned earlier, and it is now possible to produce 
quantitative non-radial oscillation models of known drift- 
ing pulsars (including PSR B0809-h74 and PSR B 0943-H0, 
Rosen fc Demorest|20lTl [Rosen fc Clemens|2008] ) , which fit 
the data at least as well as the carousel model. 

A particularly interesting feature of PSR B0809-|-74's 
driftbands is the step in subpulse phase (i.e. phase in the 
P3 direction) which i s detected at observing frequencies of 



820 MHz and a bove ( [Edwards fc S tappers| [2003| [Rosen fc 
Demorest|2011|), but is unseen at lower frequencies ([fjacker 



et al.||1975 Proszynski fc Wolszczan 1986 Edwards"^ 
Stappers| 2003). One of t he strengths of the surface oscilla- 
tion model presented in Rosen fc Demorest (2011) is that 
it can be used to explain the subpulse phase step which 
appears at high frequencies, by invoking the existence of a 
nodal line, which moves into the line-of-sight. In this paper 
we use radio observations from 14 — 5100 MHz to examine 
the driftbands of PSR B0809+74 in detail and map their 
evolution as a function of frequency, focussing in particu- 
lar on the subpulse phase step. These new observations are 
used to test current models of the pulsar magnetosphere. 

2. Observations 

Our analysis is based on newly acquired LOFAR obser- 
vations, spanning the lowest 4 octaves of the observable 
'radio window', as well as archival data from several other 
radio telescopes. LOFAR has two sets of antennas which ob- 
serve in separate frequency bands - the Low Band Antennas 
(LB As) and the High Band Antennas (HBAs). The LB As 
were used to observe the pulsar from 14 — 62 MHz. The pul- 
sar was also observed using the HBAs from 119 — 167 MHz 
and from 190 - 228 MHz. All of the LOFAR observations 
in this paper were taken using the coherent sum of the six 
stations on the 'Superterp' (central core of the array), ex- 
cept the HBA data from MJD 54896 (which were taken us- 
ing the 4-tile LOFAR test station) and MJD 55099 (which 
were taken using the incoherent sum of 3 remote stations). 
For further inf ormation on LOFAR's beamformed observ- 
ing modes see Stappers et al. (2011), and for a general 



LOFAR description see van Haarlem et al. (submitted) 

In addition to the LOFAR observations detailed above, 
we also present the results of our analysis on data at 
328 MHz and 1380 MHz f rom the Westerbork Synthesi s 



Radio Telescope (WSRT, Edwards fc Stappers 2003) 



624 MH z data from the Gia nt Meterwave Radio 'I'elescope 
(GMRT,l Gajjar et aL|20T2|), 820 MHz data fr om the Green 
Bank Telescope (GBT, Rosen fc Demorest 2011), previ- 



and new 2220 MHz data from WSRT (see [ Karupp usamyl 



et al. 



2008 for details of the system). Full details of all of 



the observations are summarised in Table [TJ 

3. Analysis 

3.1. Initial Processing 

Data from each observation were de-dispersed to a DM of 
5.75 pc cm"'^ and collapsed in frequency so that each ob- 
servation (with the exception of the one taken with the 
LOFAR LBAs) contained a single frequency channel. At 
low frequencies (below ^ 80 MHz) the pulse profile of 
PSR B0809+74 changes significantly over the wi de frac- 
tional band width of the LOFAR LBAs (see, e.g. 'Hassall[ 
et al.|2012 ), so the observation was divided into 4x12 MHz 
bands, with centre frequencies of 20 MHz, 32 MHz, 44 MHz 
and 56 MHz. Initial estimates of the pulsar's period at the 
time of ea ch observation were derived from known spin pa- 
rameters (?Hobbs et al. 2004). These estimates were op- 
timised by performing a narrow search in period space 
around the ephemeris values. Each optimised period was 
used to divide the data from the appropriate observation 
into pulse-period-sized segments and create a 'pulse stack', 
a two-dimensional array of the intensity of the source as 
a function of pulse longitude and pulse number. We used 
these pulse stacks to produce fiuctuation spectra, subpulse 
phase tracks, and P3-folded pulse stacks, in order to exam- 
ine the temporal and spectral stability of the drift patterns, 
and test models of the pulsar magnetosphere. 

3.2. Fluctuation Spectra 

To investigate the periodicity associated with drifting 
subpulses, we used the Longitud e-Resolved Fluctuation 
Spectrum (LRFS, see [Backer[[l970| . In the LRFS method, 
data are first divided into blocks of successive pulses (typ- 
ically 512 pulses per block). For each block, Fourier trans- 
forms are calculated for constant pulse longitude columns. 
The fluctuation spectra produced from all of the blocks are 
combined to produce an average fluctuation spectrum at 
each pulse longitude. If the pulsar has drifting subpulses, 
each pulse longitude bin which features a driftband will 
have a frequency associated with the distance between sub- 
sequent driftbands. These lead to 'features' in the LRFS, 
which can be used to determine P3. 



3.3. Subpulse Phase Tracks 

For most pulsars, P2 is not constant. Driftbands are often 
curved, and as a consequence, P2 varies si gnificantly with 
pulse longitude ( Edwards fc Stappers|2002 ) . Figure [l] illus- 
trates this - the black lines represent curved driftbands in a 
pulse stack, and in the highlighted driftbands, P2 is shown 
at various pulse longitudes by the grey lines. The curvature 
of the driftbands produces different values of P2 at differ- 
ent pulse longitudes. Techniques lik e the Two-dimensional 
Fluct uation Spectrum (2DFS, see Edwards fc Stappers 
2002) cannot resolve any pulse longitude-dependence, and 
ettectively provide the P2 distribution averaged over the 
whole pulse window. To investigate the shape of the drift- 
bands, we used the complex spectra produced from the 
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Table 1. Summary of the observations presented here. 



Observatory 



Centre 
Frequency (MHz) 



Bandwidth Sampling TV, 



pulses 



MJD 



Reference 



(MHz) 


iime (ms) 




48 


1.3103 


5571 


9 


1.3103 


178283 


12 


1.3103 


185628 


48 


0.6554 


33433 


38 


0.8192 


5572 


10 


0.8192 


15984 


20 


0.2288 


5696 


200 


0.1600 


478 


80 


0.4096 


13092 


120 


0.8192 


2788 


500 


0.1000 


4100 



LOFAR " 
LOFAR 
LOFAR " 
LOFAR 
LOFAR " 
WSRT 
GMRT 
GET 
WSRT 
WSRT 
Effelsberg 



38 
160 
145 
143 

209 
328 
624 
820 
1380 
2220 
4850 



55815 
54896 
55099 
55822 
56008 
51874 
55609 
54944 
52694 
56021 
51487 



This work 
This work 
This work 
This work 
This work 



Edwards & Stap 
Gajjar ot al. 


pers| (j2003| 
(2012) 


iRoscn & Dcmoi 


est (20114 


Hjdwards &: btappers ( 


2U03 



This work 



Notes. LOFAR observation IDs: L30803 L2009_11193 '^"^ L2009_14591 L30910 L53897 



LRFS. By calculating the complex phases in each pulse 
longitude bin it is possible to produce a subpulse phase 
track, which is a measure of the average shape of a drift- 
band. For a more detailed e xplanation of this process see 
Edwards fc Stappers (20021, and for a description of the 
specifi c method used in this paper see [Weltevrede et ah] 
(20T2IPI 



3.4. P3 Folding 

Subpulse phase tracks are good for probing the stability of 
driftbands from observation-to-observation, and also pro- 
vide good indicators of the shape of the driftband when the 
signal-to-noise ratio of the observations are low. However, 
in some circumstances they do not give the whole picture. 
The subpulse phase track only gives a single value for the 
subpulse phase at each pulse longitude bin, so if there are 
two overlapping driftbands present in the same pulse longi- 
tu de bin (as is suggested to be the case in PSR B0809-I-74, 
see Edwards &: Stappers|2003) ) the technique will not accu- 



rately describe the data. 

An alternative way of analysing the data, which can help 
with this problem, is P3 folding. This is done by folding the 
pulse stack (which has already been folded at the rotation 
period of the pulsar) at the P3 value derived from the LRFS, 
to produce a 2-D array containing the shape of the 'aver- 
ag e' driftband. Similar methods have been used in the past 
by lDeshpande fc RankiIIl (|200lD,| van Leeuwen et al.l (|2002| 
and Backus et al. |( |'M1| . PSR B0809-h74 is a nulling ^IHsar 
and when the pulsar goes into its null state, its drifting pat- 
tern is disrupted and the driftbands quickly go out of phase 
with each other, causing the folded driftbands to become 
smeared together. In addition, even without nulls the value 
of P3 is not stable from pulse-to-pulse and will constantly 
fluctuate around a mean value. To mitigate these effects, we 
folded the data in blocks of 256 pulses, which were added 
together using an iterative technique to maximise the corre- 
lation between subsequent blocks of subbands, by allowing 
arbitrary phase offsets between blocks. Analysis through 
the P3-folding technique is complementary to analysis with 
subpulse phase tracks, as Ps-foldcd data provides extra in- 



Note that i n this paper, we d efine subpulse phase with the 



opposite sign to Edwards & Stappers (2003) and other authors. 



We find this more intuitive, as it means that the subpulse phase 
tracks have a positive slope for subpulses which drift in the 
positive direction, so the two can be compared directly. 




Pulse Longitiide 



Fig. 1. Illustration of how curved driftbands mean that 
P2 varies with pulse longitude. The black lines represent 
the shape of simulated driftbands. The grey lines represent 
the size of P2 for the highlighted driftbands. One can see 
that P2 varies significantly depending on the gradient of 
the driftbands at a given pulse longitude. 



formation on where the power is in the pulse profile and 
can be interpreted more intuitively, but the technique does 
not work as well on low signal-to-noise data. 



4. Results 

4.1. Temporal Stability of Driftbands 



Lyne & Ashworth ( 1983 ) showed that the drift-rate of PSR 
B0809-I-74 changes regularly on very short timescales, in- 
creasing slightly after nulls, before relaxing to its pre-nuU 
speed after ~ 20 pulses. The pulsar's drift-rate has also 



been shown to v ary on timescales of several minutes by |van| 
Leeuwen et al.| ([2002), who observed a 'slow' mode which 
was stable for ~ 120 pulse^ To investigate the temporal 
stability of the drift pattern of PSR B0809+74 on longer 



^ We note that this slow mode must be quite rare and/or 
frequency dependent, as it was not observed in any of the long 
observations presented here. 



3 



Hassall et al.: Frequency-dependent Delay from the Pulsar Magnetosphere 




Fig. 2. The subpulse phase tracks of three LOFAR observations of PSR B0809+74 taken in March 2009 (white points), 
September 2009 (black points) and September 2011 (grey points). The right panel shows the same tracks with a slope of 
—25°/° subtracted and offset by ±50°, so that fine structure is easier to see. The grey line is a spline of the most recent 
observation so that data can be easily compared across epochs. The shape of the subpulse phase tracks are constant to 
within experimental uncertainties. 




Fig. 3. The subpulse phase tracks of 8x1.5 hour blocks of a LOFAR HBA observation of PSR B0809+74, showing that 
the shape subpulse phase track remains approximately constant with time (typical spread is ^ 10°). The right panel 
shows the same tracks with a slope of —25°^° subtracted, so that fine structure is easier to see. 



timescales, we used three LOFAR HBA observations span- 
ning a period of 2 years. Using the LRFS, we found P3 to 
be constant across all observations to within experimental 
uncertainties. 

The curved shape of the driftband means that P2 varies 
between 30 ms (just before the subpulse phase step) and 
43 ms (at the subpulse phase step) within a single cycle 
of subpulse phase, however the shape of the subpulse phase 
tracks remains stable on very long timescales. The subpulse 
phase tracks of the three observations are shown in the left 
panel of Figure [2] The right panel of the figure shows the 
subpulse phase tracks with a constant gradient of —25 de- 
grees (of subpulse phase) per degree (of pulse longitude) 
subtracted to show the fine structure of the subbands in 
more detail. This slope of —25°/° was chosen because it 
was the best straight-line fit to the data. From the figure, 
one can see that the subpulse phase tracks are identical (to 
within experimental uncertainties) in observations taken in 



March 2009 (white points) September 2009 (black points) 
and September 2011 (grey points). The only noticeable dif- 
ference between the tracks is their extent in pulse longitude 
- the most recent observations have more data points at the 
edge of the pulse profile compared with those taken previ- 
ously. However, this is simply a reflection of LOFAR's avail- 
able collecting area and bandwidth increasin ^ The fact 
that there are no intrinsic changes to the subpulse phase 
tracks shows that the subpulses are stable on timescales of 
months to years. 

We also used the most recent observation to investi- 
gate subpulse stability on shorter timescales. We broke the 
2011 LOFAR HBA observation into 8 pieces (each of 4096 

^ The observations from March 2009 used only 4 HBA tiles, 
the September 2009 observation used the incoherent sum of 3 
stations (each with 48 HBA tiles) and the September 2011 ob- 
servation was taken using the coherent sum of the Superterp 
(288 HBA tiles total). 
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pulses), and produced subpulse phase tracks for each piece. 
These are plotted in Figure [3j As in Figure [2] the right 
panel in the figure shows the same subpulse phase tracks 
with a constant slope of —25°/° subtracted. Again, the 
subpulse phase tracks look stable, although the spread on 
these points is noticeably larger than the spread of the long 
timescale phase tracks. We attribute the inc reased spread 
to pulse-to-pulse jitter. Helfand et al. ( 1975 1 claimed that 
the average pulse profile of PSR B0809+74 becomes stable 
after 



300 pulses, although Liu et al 
shown that the stabilisation timescales 



2012 ) have recently 



or pulsars are actu- 



ally much longer (> lO'* pulses) than expected. Driftbands 
may also need many pulses before they stabilise. This could 
explain the difference in th e size of the subpulse phase steps 
measured at 820 MHz by Rosen & Demorest (2011), who 
used relatively short integration times of 15 minutes (< 500 
pulses) in their observations. Despite the pulse-to-pulse jit- 
ter, we note that the spread in the subpulse phase tracks is 
still small - typically of order ^ 10°. Thus subpulse phase 
tracks seem to be stable on all of the timescales we have 
probed (i.e. hours, days, months and years). 



4.2. Spectral Stability of Driftbands 

We now compare the drifting behaviour of the pulsar at dif- 
ferent frequencies. From the preceding discussion, we know 
that the driftbands are stable on timescales of hours, days, 
months and years, so - although most of the observations 
we are using are not simultaneous - we can confidently as- 
sume that any differences observed with frequency are not 
caused by temporal variations. 



4.2.1. Subpulse Phase Tracks 

It has been shown previously that the shape of PSR 
B0809-|-74's driftbands vary sig nificantly as a function 
of frequency. Above 820 MHz, ([Wols zczan et al.' '1981', 

1986) [Edward s & Stappers 2003; 



Proszynski fc Wolszczan 

Rosen & Demorest 2011 1 there is a sudden step of more 



than 100" in subpulse phiase near to the peak of the pulse 
profile. But be low 500 M Hz, several authors (Edwards fc 
|Stappers|2003lat 328 MHz, [Proszynski &: Wols zczan |l986| at 



400 MHz and Backer et al.|1975 at 500 MHz) have searched 



for similar discontinuities and failed to find any evidence for 
a subpulse phase step. 

Figure |4] shows the pulse profiles and subpulse phase 
tracks obtained from the analysis of observati ons between 
14 MHz and 5100 MHz. In [Hassall et al. | j2012| , we showed 
that one of the components in PSR B0809-|-74's aver- 
age pulse profile (the component on the leading edge at 
624 MHz and below, and the trailing edge at 820 MHz and 
above) remains fixed in pulse longitude at the fiducial point, 
and the second component moves through it as a function of 
frequency. The pulse profiles in Figure |4] have been aligned 
so that this 'fiducial component' is at 0° pulse longitude. 
As in Section 4.1 the subpulse phase tracks have all had 
a gradient of —25°/° removed from them, to accentuate 
fine structure. The subpulse phase step is clearly visible at 
1380 MHz, but what is apparent is that it gradually in- 
creases in size as a function of frequency. It is visible at 
frequencies as low as 328 MHz, where it can be seen on the 
leading edge of the pulse profile. Even at 143 MHz, the sub- 
pulse phase track is asymmetric, and the leading edge looks 



like it has been disrupted. The subpulse phase step moves 
to later pulse longitudes and increases in magnitude as ob- 
serving frequency increases, reaching its maximum size of 
~ 190° at 2220 MHz. The size of the phase step is shghtly 
lower in the 4850 MHz data, although it is possible that it 
may have wrapped in subpulse phase, i.e. because the drift- 
bands recur periodically, the size of the phase step could be 
wrong by ±360°. Despite the large variation in the shape 
of the driftbands, we find (using the LFRS) that P3 is con- 
stant at all frequencies. 

4.2.2. P3 Folds 

The P3-folded pulse stacks of each of our observations are 
shown in Figure [5] The pulse stacks look radically differ- 
ent at each frequency. There are two distinct driftbands at 
the lowest frequencies, which are seemingly associated with 
the two components in the pulse profile, and the driftbands 
move closer together in both pulse longitude and subpulse 
phase with increasing frequency. This movement is most 
prominent in the LOFAR LBA observations, where the cen- 
troid of the rightmost driftband appears to arrive 10 pulses 
later at 56 MHz than it does at 20 MHz, and 23.5° earlier 
in pulse longitude. These observations (at 20, 32, 44 and 
56 MHz) were taken simultaneously, and the plots shown 
in Figure [5] are time aligned. 

The fact that one driftband remains stationary, whilst 
the other suffers a frequency-dependent delay of 10 pulses 
is contrary to the expectations of radius-to-frequency map- 
ping. Although one could de-disperse the data to a different 
DM to make it appear that the two components move away 
from each other at the same rate in the pulse profile, the 
DM of the pulsar is not large enough to account for the 
10-pulse delay in the P3 direction, which only affects one 
driftband. Thus we conclude that there are two sets of drift- 
bands - the 'fiducial driftband' (the driftband associated 
with the fiducial point of the pulse profile) appears fixed in 
both pulse longitude and subpulse phase (in strong agree- 



ment with the results ofjiassall et al.|2012|), and the other 



component suffers a frequency-dependent delay relative to 
it. 

At higher frequencies, the driftbands begin to overlap, 
and above 143 MHz it is hard to ascertain which drift- 
band is which by eye. To determine what happens above 
143 MHz, we fitted both of the folded driftbands with 2D 
gaussians using a least squares method. To alleviate some of 
the degeneracy in the region where the two driftbands over- 
lap, the shape of the fiducial driftband was held constant 
for the fit^ The half-power contours for the fitted gaus- 
sians are shown in Figure [5j The solid line represents the 
2D gaussian fitted to the fiducial driftband, and the dashed 
line represents the gaussian fitted to the other driftband. In 
the figure, the folded pulse stacks above 56 MHz have been 
aligned such that the fiducial driftband appears at a fixed 
pulse longitude and subpulse phase. The fits reveal that 
above 56 MHz second driftband continues to move through 
the pulse profile in both pulse longitude and subpulse phase. 
The moving driftband at 20 MHz arrives 30 pulses earlier 
than at 1380 MHz, whilst emission from the other driftband 
arrives at the same pulse longitude and subpulse phase at 



* We note that in the LBA observations, where the two drift- 
bands are distinct, the shape of the fiducial driftband is approx- 
imately constant. 
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Pulse Longitude (°) Pulse Longitude (°) Pulse Longitude (°) Pulse Longitude t'*) 

Fig. 5. Folded 'pulse stacks' from each of the; observations used. The solid and dashed ellipses show the half power points 
of 2D gaussians fitted to the data using a least squares method. The central observing frequency of each panel is indicated 
in the top left corner. The driftbands are repeated over 5 x P3 so that the full drift pattern is visible at all frequencies, 
yhe pulse profile (grey line), and a simulated pulse profile (black line) from each 2D fit is shown in the panel above the 
corresponding pulse stack. 
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-40 -20 20 40 60 
Pulse Longitude (°) 

Fig. 4. The subpulse phase tracks (grey points) and average 
pulse profiles (black lines) of PSR B0809+74 between 14 - 
5100 MHz. All phase tracks have had a constant slope of 
—25°/° removed. The central observing frequency of each 
panel is indicated in the bottom left corner. 
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Fig. 6. The path of the moving driftband as a function 
of frequency. The shape of the stationary ('fiducial') drift- 
band corresponding to the fiducial pulse profile component 
is shown in grey. 



all frequencies. We note that in the regions where the two 
driftbands overlap, they are not aligned. This suggests that 
the two driftbands are not part of the same system, and that 
the delay we see is not simply a result of a pulse window 
function. We attribute the fact that the driftbands appear 
to be aligned at 56 MHz to coincidence - the movement of 
the non-fiducial driftband spans 40 pulses over the fre- 
quency range observed here, so it is not unreasonable that 
the driftbands should appear to be aligned over a narrow 
frequency range. The subpulse phase step, which appears 
when the driftbands begin to overlap also suggests that the 
two driftbands are not perfectly aligned. 

Figure |6] shows the path of the moving driftband as a 
function of frequency. The black points show the position of 
the centroid of the moving driftband and the grey regions 
each represent the position and shape of the fiducial com- 
ponent's driftband. One can see that the driftband follows 
a smooth, and continuous path through the figure. Above 
1380 MHz the driftband appears to move back towards later 
pulse longitude and earlier pulse number, although it should 
be noted that this is the region in our data where the signal- 
to-noise is relatively poor and the two sets of driftbands are 
overlapping, so there is an increased degeneracy in the fits. 



5. Discussion 

In Section [4j we showed that the radio emission from PSR 
B0809+74 comes from two separate driftband structures. 
O ne driftband r e mains fixed at the fiducial point found 
in 



Hassall et al. (2012) and in subpulse phase, whilst the 
other driftband moves towards earlier pulse longitude and 
later pulse numbers with increasing frequency. The move- 
ment of the driftband is subtle above 143 MHz, and it only 
moves significantly in the LOFAR LBA band, where, un- 
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til now, narrow available bandwidths have prevented this 
effect from being seen. 

Above 143 MHz, the two driftbands overlap with each 
other in pulse longitude, and comparing the Ps-folded pulse 
stacks in Figure [5] with the position of the corresponding 
subpulse phase step in Figure |4] it is clear that the interface 
between the two different driftbands is what is producing 
the subpulse phase step. This driftband con figuration re 
places the need for the 'absorption feature' (|Bartel||1981 



Rankin|[l983! ) on the leading edge of the pulse prohle, anc 



could potentially explain the orthogonal polarisation modes 
seen in PSR B0809+74's single pulses. Certainly, the fits to 
the Pa-folded pulse stacks shown here look very similar to 
the islands o f polarisation seen in Figures 5 and 6 ofjRankm] 



et al. (2005), although unfortunately we do not yet have 
polarisation data to test this hypothesis directly. In the fol- 
lowing discussion, we attempt to reconcile our findings with 
models of the pulsar magnetosphere. 



5.1. Surface Oscillation Model 

Currently, the surface oscillation model cannot explain 
many of the phenomena shown here. One of the require- 
ments of the surface oscillation model of IClemens &: RosenI 
(2004) is that nodal lines should introduce subpulse phase 
steps of exactly 180° . The authors explained a 120° sub- 



pulse phase step seen by Edwards & Stappers (2003) in 
1380 MHz PSR B0809-h74 data by suggesting that because 
the emission we see is an average over a finite frequency 
range, the phase step is "washed out". This is a valid ar- 
gument for why one frequency would have a phase step not 
exactly equal to 180°, but it is unlikely that the subpulse 
phase step would increase so systematically with frequency 
(see Figure l4| , despite the differing bandwidths in each ob- 
servation, ifthis were the case for all frequencies. 

A second requirement of the model is that the spacing 
of the components should follow the same distribution as 
a spherical harmonic sampled along a single line-of-sight. 
We can rule this out qualitatively for PSR B0809-I-74 by 
noting that the size of the fiducial component is approxi- 
mately constant at all frequencies, and the spacing between 
the components in pulse longitude varies (smoothly) from 
^ — 5° to ^ -1-35°. The model also cannot explain how the 
subpulse phase step begins on the leading edge of the pulse 
profile at low frequencies and moves through the central 
component, appearing on the trailing edge of the compo- 
nent at high frequencies. 

5.2. Carousel Model 



but to account for what we observe here, absorption needs 
to selectively effect only the leading component at low fre- 
quencies, and the trailing component at high frequencies, 
which seems imlikely. Even using a distorted polar cap (as 
discussed in Arendt & Eilek 1998), it is impossible to re- 
produce the pulsar's profileevolution . In the latest E x B 
models (e.g. Timokhin fc Arons|2012 1, the currents and po- 
tential differences in the magnetosphere can be asymmetric 
along the line of sight, van Leeuwen & Timokhin (2012) 
used this fact to explain the two different drift speeds seen 
in PSR B0826-34. Despite this, such models do not change 
the geometry of the system, and therefore cannot explain 
the two sets of driftbands, or the movement of one driftband 
through the other as a function of frequency. 

Even invoking two carousels (one for each driftband), 
the emission from one driftband must undergo a strongly 
frequency-dependent delay in subpulse phase whilst all of 
the emission from the other must arrive almost simulta- 
neousljj^ If we take this delay at face value, as a physical 
30-pulse delay, the size of the region needed to produce such 
a delay by light-travel time is 10^ km, or ~ 180 times the 
radius of the light cylinder. Alternatively, we can interpret 
the changing drift pattern as being due to the one pulse- 
profile component moving to later pulse longitudes, being 
modulated by an underlying driftband. This scenario still 
requires a delay of 40° in pulse longitude (~ 0.14 seconds), 
corresponding to a light-travel time of 43000 km, roughly 
70% of the radius of the light cylinder. Both interpretation s 
are in stark contrast with the results of Hassall et al. ( 2012 ) , 
where we showed that all of the emission from the fiducial 
component in the pulse profile between 46 MHz and 8 GHz 
must originate from within a region smaller than 400 km. 



5.3. Birefringence 



In Hassall et al. (2012), we suggested that birefringence 
could be responsible for the evolution of the pulse profile 
with frequency, but the fact that the shape and pitch an- 
gle of the moving driftband changes with frequency sug- 
gests that it is unlikely to be caused by a simple, axis- 
symmetric propagation effect. If it were, the driftbands 
would move relative to each other but the shape and pitch 
angle of the driftbands would remain constant at all fre- 
quencies. Perhaps refraction could account for the observed 
frequency-dependence when coupled with an asymmetric 
plasma distribution in the magnetosphere, although the 
frequency-dependent delay also seems large for a propa- 
gation effect which must occur in a region smaller than the 
light cylinder (10^ km). 



The behaviour seen here also cannot currently be explained 
using the carousel model. Using a carousel, it is impossible 
to recreate the two components which are distinct below 
143 MHz and overlap at higher frequencies. The carousel 
model has an inherent symmetry to it, and when compo- 
nents in the pulse profile move apart they should do so 
symmetrically to fit with the hollow cone described by the 
inner edges of the dipolar magnetic field. The average pulse 
profile of PSR B0809-I-74 (particularly its spectral evolu- 
tion) is very asymmetric. Most scenarios which are able to 
explain the movement of the component require one half of 
the carousel to be missing. Abs orption has been used in th e 
past to explain this asymmetry ( Bartel|1981 Rankin 1983 1, 



6. Conclusion 

We have shown that both the shape, and the separation 
between the driftbands (Fg) of PSR B0809+74 are stable 
on timescales of '-^ years. But, whilst P3 does not change, 
the shape of the driftbands is extremely variable with fre- 
quency. Emission from PSR B0809-I-74 is composed of two 
separate driftbands. The emission from one of the drift- 
bands arrives at the same pulse longitude in t he same pulse 
at all frequencies (the fiducial point found in Hassall et al. 



^ Note that this argument is equally valid for rejecting a model 
with two surface-oscillation modes. 
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2012 1, whilst emission from the other driftband changes lo- 



cation in pulse longitude and is delayed by tens of pulses 
in subpulse phase at high frequencies. The subpulse phase 
step, which is visible at high frequencies, is attri buted to the 
interface of the two driftbands (as suggested by |Edwards fc 
Stappers||2003[ ). 

The carousel and the surface-oscillation models are not 
emission mechanisms per se, and so, should not be expected 
to explain frequency-dependent effects. However, the fea- 
tures seen in the spectral evolution of the driftbands shown 
in our observations seem incompatible with current incar- 
nations of either model. We are, as of yet, unable to find 
any emission mechanism or geometry to explain these fea- 
tures satisfactorily. Similar studies on the spectral evolu- 
tion of subpulse phase in other pulsars at low frequencies 
will be done to determine whether PSR B0809-I-74 is a rare 
case, or if this type of frequency evolution is common to 
the pulsar population. More data could also provide fur- 
ther clues to the mechanism causing the effect. In partic- 
ular, low-frequency polarisation data would be very useful 
in determining whether this phenomenon is related to the 
orthog onal polarisation modes, as seen by [Rankin et al. 
(2005) and others. It could also be used to extract more 
information about the conditions of the magnetosphere in 
the regions where each of the driftbands originated. 
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